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Near-threshold rj production in pp collisions 

Qi-Fang Lii and De-Min lQ 

Department of Physics, Zhengzhou University, Zhengzhou, Henan 450001, China 

We study near-threshold meson production in pp collisions within an effective 
Lagrangian approach combined with the isobar model, by allowing for the various 
intermediate nucleon resonances due to the ^ r, p, and p-meson exchanges. It is 
shown that the p-meson exchange is the dominant excitation mechanism for these 
resonances, and the contribution from the IV*(1720) is dominant. The total cross 
section data can be reasonably reproduced, and the anisotropic angular distributions 
of the emitted p meson are consistent with experimental measurements. Besides, the 
invariant mass spectra of pp and pp explain the data well at excess energy of 15 MeV, 
and are basically consistent with the data at excess energy of 40 MeV. However, 
our model calculations cannot reasonably account for the two-peak structure in the 
pp distribution at excess energies of 57 and 72 MeV, which suggests that a more 
complicated mechanism is needed at higher energy region. 

PACS numbers: 13.75.Cs; 14.20.Gk; 13.30.Eg 


I. INTRODUCTION 

The properties of hadrons are subject to the behavior of quantum chromodynamics (QCD) 
in the non-perturbative region, so the study of hadrons is important to deepen the under¬ 
standing of the non-perturbative properties of QCD. Meson production reactions in nucleon- 
nucleon collisions near threshold are a good platform to obtain new information on the hadrons 



p meson production in NN collisions is generally assumed to occur predominantly through 
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re-scattering of the intermediate nucleon resonances caused by the meson exchanges. For this 


basic mechanism, it is not yet c 


ear which of the possible nucleon resonances plays the dominant 


role. For example, in Refs. 15-|l9l, 22-28] it is suggested that the iV*(1535) is dominant, while 
in Ref. 33] it is found that the A r *(1520) is dominant and the iV*(1535) contribution is small 


due to the strong destructive interference among the exchanged mesons. Also, it is not yet 


clear which o 
Refs. 


15 


221 , 


the possible meson exchanges plays the most important role. For example, in 
32| it_is claimed that the pseudoscalar mesons 7 r and 77 exchanges are dominant, 


whereas in Refs 
Ref. 


3,Q, 


24 


25] it is suggested that vector meson p exchange is dominant. In 


27] it is found that both the 7 r exchange and the p exchange can describe the cross sections 


well. The study on the !V*(1535)lVp coupling in the framework of an effective Lagrangian 
approach manifests that the value of gN*(i 535 )Nf> is strong [35|, which favors the importance of 
the p meson exchange in NN collisions. 

In Ref. [34|, the final state interaction (FSI) enhancement factor is considered and it is found 
that the measured pp and pp effective mass spectra can be well reproduced by allowing for a 
linear energy dependence in the leading 3 P 0 —> 1 S 0 ,s partial wave amplitude. 

In Ref. [ 8 ] it is suggested that the higher partial waves may be important even at 15.5 MeV. 
Besides the A r *(1535), the p meson may also couple strongly to other higher resonances. The 
large branching ratio and the small phase space for the IV* (1720) —> Np also suggest that the 
N* (1720) Np coupling is strong. 

With the inspiration of the factors mentioned above, we shall restudy the pp —> ppp reaction 
in an effective Lagrangian approach combined with the isobar model. The combination of the 
effective Lagrangian approach and the isobar model turns out to be a good method to study 


hadron resonance production in the nN, NN, and KN scattering 22 


32 


33 


35 


40]. In the 


present work, we assume that the near threshold 77 meson production in proton-proton collisions 
is through the intermediate IV* (1535), IV* (1650), IV* (1710), IV* (1720), and the nucleon pole 
caused by the i r, 77, and p mesons exchanges. The proton-proton FSI and proton -77 FSI are also 
considered. 

This work is organized as follows. The basic formalism and ingredients used in our model are 
given in Section II. The numerical results and discussions are given in Section III. A summary 
is given in Section IV. 
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II. FORMALISM AND INGREDIENTS 


The basic tree level Feynman diagrams for the pp ppr) reaction are shown in Fig. [H 



FIG. 1: Feynman diagrams for the pp —» ppp reaction. 


37j]: 


The interaction Lagrangians for the nNN, r]NN and pNN couplings can be written as 


36 


£~kNN — — ig-KNN'^Nlb'T ■ ft'll’N, 
£ijNN = ^igrtNN^Nlb^N-, 




K 


£ P NN = dpNN'IpN ( 7/u + ^ ' P^^N- 


( 1 ) 

( 2 ) 

(3) 


The effects of the non-point-like structures of exchanged mesons are taken into account by 


introducing the following off-shell form factors in the A4NN vertexes 

F nn (k 2 ) — 

r M K^MJ ~ 


41 


431 


AL — m 


M 

A 2 


k 2 


(4) 


where M denotes the exchanged meson; n — 1 for tt°NN and r/NN vertexes, n — 2 for the 
p°NN vertex; Izm, m-M, and A m are the 4-momentum, mass, and cut-off parameter, respectively, 


in our calculations are: q 2 y V /47T 
-481, Ap = !-85 GeV 36„ [33], and 


for the exchanged-meson M. The relevant parameters usee 
= 14-4, gf lNN / 47T = 0.4, g 2 p nn /Att = 0.9, k = 6.1 3fJ, \37 
A n = Ap = 0.8 GeV. 

The following interaction Lagrangians involving the nucleon resonances N* can be obtained 
within a Lorentz covariant orbital-spin (L-S) scheme for the N*NM couplings 42]: 


AriVJV*(1535) — *#7riV7V*(1535)f/ ; iV'r • 7T0JV*(1535) + h.C., 


(5) 








































£-t 1 NN*( 1535) — WrjNN* (1535) $NV$N* (1535) + h.C., 
f-'pNN* (1535) = WpNN* (1535)$ Nl5 (ip — ^ ' P^ (Pp)$N* (1535) + h.C., 

£ttAT iV*(1650) = ^9txNN*(1050)$N? ' 7I"0iV*(165O) + h.C., 
A/AW*(1650) = ^9r]NN* (1650) $NV$N* (1650) + h.C., 

^pAfiV*(1650) = WpNN*(1650)$Nl5 (ip ~ ? ' P^ (P p)$N* (1650) + h.C., 




r 


^wNN* (1710) — WttNN*(1710)$N15? ' 7T0iV*(171O) + ^- c -) 

-^-'j/JVTV* (1710) = —^9pN N* (1710)$ Nl59$ N* (1710) + h.C., 

FpNN* (1710) — — 9pNN*(1710)$N (ip + T^ - &T ■ P^$N* (1710) + h.C., 
h-'trNN*(1720) = 9irNN* (1720)$N? ’ d tl 'K$N*(1720)p + h.C., 

FpNN* (1720) = g V NN* (1720)$ N 9^ 9$ N* (1720)/u + h.C., 

£-pNN*(1720) ~ 9pNN*(1720)$N15? ‘ P^$N*(1720)p + h.C.. 

At the N*NM vertexes, the monopole form factors are employed: 


N (k 2 1 - 
"m — 


A 


*2 

M 


m 


M 


A *2 _ 7,2 • 

iv M 


(6) 

(7) 

( 8 ) 
(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 


With the effective Lagrangians listed above, the partial width of the nucleon resonance W* 
can be derived. From the experimental data on the partial width of the corresponding nucleon 
resonance, the N*NM coupling constants can be obtained. For the W* resonance below the Wp 
threshold, the partial decay width F N*->Np->Nirn is employed to determine the N*Np coupling 
constant (361.1371], The relevant coupling constants and cut-off parameters are listed in Table H 



A>r 

51- 

-53] 


F N (q 2 ) = 


A^r 


F N $q 2 ) = 


with A/v = 1.0 GeV and A N * = 2.0 GeV 


a n + (v 2 - m n ) 2 ’ 

_A AT*_ 

A f'. + tf-M*.)*’ 


(18) 

(19) 


The meson propagators used in our calculation are: 
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TABLE I: Relevant parameters of the nucleon resonances used in our calculation. The widths and 
branching ratios are taken from the PDG [5C|. 


Resonance Width(GeV) Decay channel Branching ratios g 1 j 47r Cut-off(GeV) 


A*(1535) 

0.15 

Nn 

0.45 

0.037 

0.8 



Nq 

0.42 

0.28 

0.8 



Np 

0.02 

5.55 

0.8 

A*(1650) 

0.15 

Nn 

0.70 

0.052 

1.5 



Np 

0.10 

0.036 

1.5 



Np 

0.01 

0.0064 

1.5 

2V*(1710) 

0.1 

Nn 

0.125 

0.072 

1.5 



Nq 

0.20 

0.97 

1.5 



Np 

0.15 

0.019 

1.5 

1V*(1720) 

0.25 

Nn 

0.11 

0.11 

1.5 



Nq 

0.04 

0.35 

1.5 



Np 

0.775 

635.11 

1.5 


G?%) = ~i 


. (- KK/ k l 


kp — m? p 


The propagators of the N* resonances can be written as 

4 + m n * 


for spin- 7 ; resonances, and 


o N *(q ) = 


os', fa) = 


(p — + iM N *T jy* 

_ ~Py.v{g) _ 

q 2 — M'n* “I" iM N *Y N * 


with 


Phv(q) — — (4 P Mn*) 


1 

3 


9^ r/lyTiv 


3 Mi 


N 1 


\ln<h> - 7 


for spin -1 resonances. 

As usual, in Eq. (1 22 ]) t 
IV (1535) (s), is employed 


( 21 ) 


( 22 ) 


(23) 


(24) 


re 


54 


bllowing energy-dependent total width of the A^*(1535) resonance, 


T 


7V*(1535) 


» = r 


Pttn{ s ) r PriN{ s ) 

N* (1535)—>W7r-7TT2-T + 1 7V*(1535)->AV 

PttN { M n *( 1535W 


Pr)N(M AT* (1535)) 


(25) 
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where the two-body phase space factor p 7r ( J? )Ar(s) i 


is 




2 # 


(r))N ' 


(s) _ y/[s - (m N + m^ ( , ? )) 2 ][s - (m N - m,;,)) 


(26) 


y/s s 

The invariant amplitude for the nucleon pole N or nucleon resonances N* can be expressed 


as 


m n / n *= m 


N/N* 


i=rr, rj, p 


M 


N/N* 


^ M l 

j=a, b, c, d, 


N/N* 

i 


(27) 


(28) 


where p a = pd = 1 and rjb = p c = —1. The explicit expressions of can be derived 

straightforwardly according to the Feynman rules. For example, A in,a 153 ^ can be written as 


M-nJ 1535 ' > — gTTNN9TrNN*(1535)9riNN*(1535)F^ N (kl)F^ l 1535 )^ (fc^)F/v*(i535) (<? 2 

x G v (k v )u{p 3 , s 3 )-/ 5 u(p 2 , s 2 )u(p 4 , s 4 )G N * { 1535 ) (q)u(p u si), 


(29) 


where s* (i = 1,2, 3,4) and pt (i = 1,2, 3,4) represent the spin projection and 4-momentum of 
the two initial and two final protons, respectively. 

The pp FSI and pp FSI are taken into account by introducing the following enhancement 
factors 

kpp + 


Fpp(kp P ) 


kpp — ia ’ 


(30) 


Fpp ( kpp) 


1 — ik pv a ’ 


(31) 
y. The 

relevant parameters are: a = 0.1 fm _1 , j3 = 0.5 fm _1 55], and a = (0.487 + zO.171) fm 52,153]. 


where k pp and k prj are the internal momenta of the pp and pp subsystems, respective 


The overall final state interaction is the product of these enhancements 

F fs i Fp P ( kpp) Fp P ( k PP3 ) Fpp ( k VPi ), 


29 


55|: 


(32) 


where p 3 and p 4 denote the two final protons. 

With the modular square of the full invariant amplitude |-M| 2 = Y A4 A 7- v * p FSJ | 2 , the 

N,N* 

differential cross section for the pp —> ppp reaction can be written as 

7 / . n 1 m P i A A \ 2 m p ( Pp 3 rn p d 3 p 4 d 3 p 5 1 4 , . , 

dcr(pp ^ ppp) =-—2^2^\M\ —---- ttft(Pi +P 2 ~P3 ~Pa ~P 5 ), (33) 


F.> 


Ea 2Ek 2 
























where F is the flux factor 
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F=(27r) 5 v /(p 1 -p 2 ) 2 -m4. (34) 

The factor | before the 5 function in Eq. (fool) comes from the two identical protons in the final 
states. The interference terms between different resonances are ignored. 


III. NUMERICAL RESULTS AND DISCUSSIONS 


With a Monte Carlo multi-particle phase space integration program, the total cross section 
versus excess energy £ up to 80 MeV, the invariant mass spectra, angular distributions, and 
Dalitz plots at excess energies £ = 15, 40, 57, and 72 MeV for the pp —> ppp reaction are 
calculated. 

The total cross section is shown in Fig. [2] together with the experimental data. Our results 
agree fairly well with the experimental data. From Fig. [21 one can see that the contributions 
from the t-channel p and 7T-meson exchanges are important and the p exchange plays the 
dominant role, but the contribution from the ? 7 -meson exchange is negligible. Fig. [2] also shows 
that the contributions from the IV* (1720) and A* (1535) are important and the A* (1720) plays 
the dominant role. The contribution of the A* (1535) is smaller than that of the A* (1720) due 
to the strong destructive interference among the exchange mesons, which is similar to the result 
]. The contributions from the A* (1650) and A* (1710) are negligible. 


of Ref. 


The invariant mass spectra, angular distribution, and Dalitz plot at excess energy £ = 15 
MeV are shown in Fig. [3] together with experimental data. The measured pp and pp invariant 
mass spectra and the angular distribution of p are reproduced well. From Fig. [3] (a) and (d), 
one can see that the pp FSI plays an important role. 

The invariant mass spectra, angular distribution, and Dalitz plot at excess energy £ = 40 
MeV as well as the experimental data are shown in Fig. [4] For the invariant mass spectra of 
proton-proton and proton-?/, the theoretical results are in agreement with the experimental data 
except for those near the proton-proton (proton-?/) threshold. This small discrepancy indicates 
the pp FSI used in our calculation may be somewhat strong in this region. 

For the angular distribution of the emitted p meson in the overall c.m. frame, there are two 
groups of data which do not agree with each other |s, 57]. One is isotropic 57], while the other 
is anisotropic js], as shown in Fig. [4] (c). Our result indicates that the angular distribution 










FIG. 2: Total cross section vs excess energy e for the pp_-^r ppp reaction from the present calculation 
(solid curves) are compared with experimental data [243, 15614591]. (a): The dashed, dotted, and 


dashed-dotted lines stand for contributions from the ^ r, p, and ?/-meson exchanges, respectively, (b): 
The dashed, dotted, short-dotted, dashed-dotted, and dot-short-dashed lines stand for contributions 
from the N, IV*(1535), IV*(1650), 1V*(1710), and IV*(1720), respectively. 

of the 7? meson is anisotropic, consistent with the data from Ref. 8]. As pointed out by 
Ref. js, 17, 58], the anisotropy is probably due to a mainly destructive interference between 
the dominant p exchange and 7r exchange. It is interesting to point out that the IV* (1535) 


dominant interpretations 


22 


28 


60] give almost isotropic angular distribution of the 77 at this 


region except that Ref. 27] gives the anisotropic angular distribution of the 7/ by allowing for 


contributions from baryonic and mesonic currents. 

The invariant mass spectra, angular distribution, and Dalitz plot at excess energy e = 72 
MeV as well as the experimental data are shown in Fig. [5] The experimental data shown in 
Fig. [5] (a) indicate the pp FSI should be rather weak, so the pp FSI is ignored in this energy 
region. This rough procedure has been used in double-pion production in nucleon-nucleon 


6l| . Our pp invariant mass 


collisions and the results turn out to be considerably improved 
spectrum can reasonably account for the data. 

The two-peak structure in the proton-// distribution cannot be reproduced in our calcula¬ 


tion, which is similar to the result from Ref. 


33]. This suggests that the structure in the pp 


distribution cannot be simply interpreted by the N* (1535), N* (1650), N* (1710), and N* (1720) 
resonances, and a more complicated mechanism is strongly called for. 
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FIG. 3: Differential cross sections (solid lines) and Dalitz plot for the pp —>• ppp reaction at the excess 
energy of e = 15 MeV compared with the experimental data 56j, |57] and phase space distribution 
(dashed lines), (a) Distribution of the square of proton-proton invariant mass; (b) Distribution of the 
square of proton-?] invariant mass; (c) Angular distribution of the emitted ?y meson in the c.m frame 
of the total system; (d) Dalitz plot. 

Our angular distribution of the rj at £ = 72 MeV again indicates that the r] distribution 
is anisotropic, consistent with the data from Ref. 8]. To our knowledge, there is as yet no 
theoretical paper for addressing the angular distribution of the ij at this region. 

The invariant mass spectra, angular distribution, and Dalitz plot at excess energy e = 57 
MeV as well as preliminary experimental data are shown in Fig. [6l Similar to the case at excess 


energy £ = 72 MeV, the preliminary data of Ref. 


62] show that the two-peak structure appears 


in the proton-// distribution and the angular distribution of the ?/ is anisotropic. Our predicted 
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FIG. 4: Same as Fig. [3] 
(squares) and Ref. [8] 


3] but at excess energy of e = 40 MeV. Experimental data are taken from Ref. 57] 
(dots). 


angular distribution for the ij at e = 57 MeV is anisotropic, consistent with the data. However, 
the pi] invariant mass distribution shows large differences between the present model and the 
experimental data. 

It is noted that our present model does not include the higher partial waves for the pi] 
FSI. As pointed out by Refs, js, 63], the experimental data indicate the higher partial waves at 
higher reaction energies could be important. The contributions of the higher partial waves being 
neglected in the current model calculations may cause the discrepancy between the experimental 
data and our model of the proton-7/ distribution at excess energies e = 57 and 72 MeV. 
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FIG. 5: Same as Fig. [3] but at excess energy of e = 72 MeV, and with pp FSI ignored. Experimental 
data are taken from Ref. 


IV. SUMMARY AND CONCLUSIONS 

In this paper we have calculated the pp —» ppij reaction within an effective Lagrangian 
approach combined with the isobar model. Our model calculations can reasonably reproduce 
the total cross sections up to excess energy of 80 MeV. 

It is shown that for the pp —> ppr] reaction, the contribution of the p-meson exchange is 
larger than that of the 7r-meson exchange, and the contribution of the N* (1720) is larger than 
that of the N* (1535). 

Also, the same cut-off parameters for the Af*(1650), IV* (1710), and V*(1720) resonances 
are used, which makes it suitable to investigate the relative contributions of the IV* (1650), 
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FIG. 6: Same as Fig. [3] but at excess energy of e = 57 MeV, and with pp FSI ignored. Experimental 
data are taken from Ref. [621 ]. 


IV* (1710), and N* (1720) resonances. Our results show that the contributions from the N* (1650) 
and N* (1710) are negligible. 

Our calculations can reasonably explain the measured pp and pp invariant mass spectra at 
excess energies s = 15 and 40 MeV, but fail to explain the two-peak structure in the proton-7/ 
distribution at excess energies e = 57 and 72 MeV, which suggests that in the higher energy 
region, a more complicated mechanism is needed. 

We give the anisotropic angular distribution of the r\ at £ — 40, 57 and 72 MeV, consistent 


with the data from Refs. 


62]. This favors the interpretation that the interference between 


the p exchange and n exchange is mainly destructive. 
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